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1 .  INTRODUCTION 


A  knowledge  of  the  fluid  dynamic  and  optical  properties  of  explosion- 
generated  dust  clouds  is  a  prerequisite  to  the  realistic  assessment  of  electro- 
optical  (EO)  system  battlefield  performance  degradation.  Aerodyne  Research, 
Inc.  has  constructed  a  model  --  DIRTRAN  (Disturbed  Infrared  Transmission)  — 
which  describes  those  time/space  dependent  cloud  properties  as  a  function  of 
explosive  charge,  soil  characteristics,  atmospheric  conditions,  and  EO  system 
operating  wavelengths.  This  report  describes  the  structure  and  features  of 
the  components  of  the  model,  presents  selected  results  of  its  validation 
against  photographic  and  optical  data  obtained  In  the  DIRT  I  and  Graf  II- 
Winter  field  trials,  and  discusses  issues  which  must  be  resolved  by  future 
theoretical  and  experimental  studies. 

The  components  of  the  model  are  discussed  in  sequence,  namely:  (1)  the 
dust  source  term  and  initial  cloud  formation;  (2)  the  diffusion  of  the  initial 
cloud  under  the  action  of  buoyant -rise/wind-driven  diffusion  and  gravitational 
settling  of  the  dust  particles;  (3)  the  time  dependent  extinction  of  visible, 
infrared,  and  millimeter  wavelength  radiation  by  the  dust  cloud  along  pre¬ 
scribed  lines  of  sight. 

A  computer  code  —  DIRTRAN-I  —  which  embodies  the  DIRTRAN  model  and  is 
compatible  with  the  Electro  Optics  System  Atmospheric  Effects  Library 
(EO-SAEL) ,  is  described  in  Ref.  1  and  has  been  delivered  to  the  Army  Atmo¬ 
spheric  Sciences  Laboratory. 
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2.  CHARACTERISTICS  OF  DUST  SOURCE  TERM  AND 
INITIAL  DUST  CLOUDS 


2 . 1  Airborne  Dust  Source 

The  quantitative  determination  of  the  aerosol  source  term  depends  on 
two  major  steps:  (1)  the  mass  of  aerosol  injected  into  the  atmosphere; 

(2)  the  distribution  function  of  that  mass  as  a  function  of  dust  particle 
size.  The  practical  significance  of  accurate  analyses  is  illustrated  in 
Figs.  1  and  2. 

Figure  1  shows  the  results  of  transmission  tests  in  the  8-12  pm  band 
performed  by  the  Army  Night  Vision  &  Electro-Optics  Laboratory  (NV&EOL) , 

Ft.  Belvoir,  VA.  This  test  —  known  as  TAP ATS  (Threat  Artillery  Preparation 
Against  Tank  Sights)  was  conducted  in  July  1978  at  Ft.  Knox,  Ky.  The  Trans¬ 
mission  data  shown  are  the  results  of  viewing  through  a  155  mm  HE  barrage  of 
intensity  0.7  rounds/km/sec.  Notable  in  the  figure  is  the  fact  that  from  the 
first-round  impact  to  several  tens  of  seconds  after  the  barrage  ended,  the 
transmission  was  essentially  zero  (at  best  4%  at  regular  intervals) .  In 
Fig.  2  is  plotted  the  transmission  during  a  barrage  conducted  at  Grafenwohr, 
Germany  in  November  1978  by  NV&EOL.  In  this  case,  the  barrage  had  an  intensity 
of  1.5  rounds/km/sec.  Here,  attenuation  due  to  background  fog  (as  determined 
prior  to  the  barrage)  has  been  removed  from  the  transmission  data  so  that  what 
is  shown  in  the  figure  is  strictly  the  loss  in  transmission  to  the  barrage. 

It  is  evident  from  the  figure  that  the  transmission  was  considerably  in  excess 
of  20  to  40%  for  several  tens  of  seconds  at  several  different  times  during  the 
barrage.  We  point  out  that  relative  to  the  Graf  II  case,  the  TAPATS  test  had 

a.  Comparable  soil  moisture  content 

b.  Half  the  barrage  intensity 

c.  50%  greater  windspeed 

d.  Greater  atmospheric  diffusivity  (approximately  Pasquill  Category  B 
versus  C  for  Graf  II). 
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Transmission  in  the  8-12  pm  Band  During  a  155  mm  HE  Barrage  at 
Ft.  Knox,  Ky.  For  this  barrage  intensity  and  atmospheric  con¬ 
dition,  essentially  zero  transmission  (at  most  U7,  at  regular 
intervals)  was  observed  for  the  duration  of  the  barrage  and 
several  tens  of  seconds  after  it.  (Data  from  the  TAFATS  test 
conducted  by  the  Night  Vision  f.  Electro-Optics  Laboru'orv, 

Ft.  Belvoir,  VA,  July  1978.) 
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Figure  2.  Transmission  in  the  8-12  um  and  0. 8-1.1  pm  Band  During  a  155  me  KE 
Barrage  at  Grafenwohr,  Germany.  In  this  situation,  transmission 
in  excess  of  20  -  40%  occurred  at  numerous  times  during  the 
barrage.  (Data  from  the  Graf  II  test  conducted  by  the  Night 
Vision  &  Electro-Optics  Laboratory,  Ft.  Belvoir,  VA,  November  1978.) 


y  *!**.*> 


? 


Despite  this  reduced  barrage  intensity  and  increased  aerosol  dispersion,  the 
TAP ATS  test  still  had  an  order  of  magnitude  less  transmission.  The  reason 
for  the  enhanced  transmission  in  Fig.  2  vis-a-vis  Fig.  1  will  become  manifest 
upon  the  following  discussion  of  the  two  basic  steps  in  the  determination  of 
the  dust  source  term. 


2.1.1  Mass  of  Airborne  Ejecta 

There  is  an  extensive  background  (see  for  example  Ref.  2)  of  theoretical 

analysis  and  empirical  correlations  for  cratering  mechanics  and  phenomenology. 

(3) 

Dimensional  analysis  readily  shows  that  the  governing  parameters  are  the 
density  of  the  soil,  its  elastic  and  plastic  moduli,  the  gravity  acceleration, 
and  the  energy  coupled  into  the  ground.  In  fact,  the  crater  volume  can  be 
calculated  from  numerical  integrations  of  the  conservation  equations  for  mass, 

(4) 

momentum,  and  energy  together  with  the  equation  of  state,  to  obtain  results 

in  reasonable  agreement  with  observations.  However,  the  execution  of  these 

calculations  is  quite  burdensome  and  the  number  of  available  examples  is  limited. 

Therefore,  one  cannot  extract  the  systematic  dependence  of  the  crater  volume  on 

elastic  and  plastic  properties  of  various  soil  types.  The  difficulty  has  beer. 

remedied  in  practice  by  the  use  of  empirical  correlations  of  the  very  extensive 

crater  volume  data  basis^  In  these  correlations,  the  scaled  crater  radius 

r^  and  depth  d^  are  predicted  as  a  function  of  a  scaled  charge  depth  X^  and 

qualitative  soil  characterizations  (e.g.,  dry  and  moist  cohesive).  Analytical 

fits  to  the  dependence  of  crater  dimensions  r  and  d  as  a  function  of  X  and 

c  c  c 

soil  category  have  been  presented  in  Ref.  6  and  are  included  in  the  PIRTRAN-I 
code.  Specifically,  for  bare-charge  explosions 


an  +  a  X 

0  1  c 


+  a„X2  +  a,^ 
2  c  3  c 


d 

c 


b.  +  b  x  +  b.x2  +  b,x3 
0  1  c  2  c  3  c 


+  b.X 


4 
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where 


r  ■  (r  /W1^3)  -  scaled  crater  radius  in  units  of  (m/kg1^3) 

c  c 

dc  “  (d^/W^*3)  ■  scaled  crater  depth  in  (m/kg^*3) 

Xc  *  (dfa/W^3)  •  scaled  explosion  depth  in  (m/kg^3) 

W  -  energy  release  of  charge  in  equivalent  kilograms  of  TNT 

These  coefficients,  Sq,  ...,  and  b^,  ....  arelisted  in  Tables2-1  and2-2. 
Their  average  values  for  the  various  soil  categories  recognized  in  Ref.  5  are 
shown  in  Figs.  3  and  4.  The  usefulness  of  these  curve  fits  is  apparent  from 
the  mean  crater  volumes  observed  in  the  DIRT  I  series  (Fig.  5)^  They  cor¬ 
respond  to  a  soil  intermediate  between  "dry  sandy  soil"  and  "dry-to-moist 

sand."  Since  this  site  is  characterized  as  "brown  silty  sand"  with  an 

(8  ) 

average  of  7 %  moisture,  the  qualitative  description  is  reasonable.  Similarly, 

for  the  Graf  II-Winter  test,  the  mean  crater  volume  corresponds  to  a  "wet  sand" 

type  soil,  which  appears  consistent  with  a  description  of  "  loamy  sand/sandy 

(9) 

loam"  with  an  average  moisture  content  of  20%,  The  coefficients  appropriate 
to  these  soils  are  incorporated  in  the  DIRTRAN-I  code  as  descriptive  of 
"Desert"  and  "European"  conditions,  respectively. 

The  tabulated  coefficients  must  be  unique  functions  of  the  elastic/plastic 
parameters  indicated  by  dimensional  analysis.  Thus,  it  would  seem  highly 
desirable  to  replace  the  qualitative  characterization  of  soils  with  a  statement 
of  mean  values  of  elastic/plastic  parameters  appropriate  to  each  category. 

Static  laboratory  measurements  of  an  adequate  number  of  soil  samples  for  each 
category  can  yield  the  values  of  such  parameters.  Thereafter,  the  determina¬ 
tion  of  crater  volume  would  be  uniquely  related  to  the  measurement  of  specific 
elastic/plastic  properties.  It  should  be  noted  that  crater  measurements  at  a 
single  test  site  (e.g.,  see  Refs. 7-9  for  the  DIRT  I  and  Graf  II-Winter  tests) 
show  scatter  in  crater  volume  data  of  typically  +20%  (cf.  Fig.  5).  It  would 
be  interesting  to  correlate  this  scatter  with  that  exhibited  by  corresponding 
measurements  of  the  elastic/plastic  moduli. 
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Table  2-1.  Coefficients  for  Average  Scaled  Crater  Radius 


^^s^oef  ficient 

Category 

ao 

al 

a2 

a3 

I.  II 

0.271 

-0  684 

0.390 

0.886 

III 

0.386 

-  0.849 

0.367 

0.993 

IV 

0.503 

-0.954 

0.450 

1.19 

V,  VI 

0.629 

-1.08 

0.264 

1.12 

VII 

0.806 

-1.28 

-0.178 

i 

0.852 

j 


Table  2-2.  Coefficients  for  Average  Scaled  Crater  Depth 


^  ^Coefficient  1 

b  . 

I 

0 

Category 

-0.477 

-0.571 

-0.840 

-1.17 

-1.49 

-1.82 


Ca  tegory 

Descript  ion 

I 

Rock 

II 

Dry  Cohesive  Soils  (Lightly  Cemented’' 

Ill 

Dry  Sandv  Soils 

IV 

Drv-to-Moist  Sand  and  Frozen  Ground 

V 

Soft  Rock 

VI 

Wet  Sand,  Moist  Cohesive  Soils,  and.  Ice 

VI  I 

Wet  Cohesive  Soils  (Not  Saturated!  and  Sr.ov 

Average 


Scaled  Burst  Depth  (m-kg-^) 


Figure  3.  Average  Scaled  Crater  Radius  Versus  Scaled  Burst  Depth  for 

Different  Soil  Characterizations  (Ref.  5)  Using  Coefficients 
From  Tables  I  and  II  (Ref.  6) 
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Average  Scaled  Crater  Depth  (m-kg 


Scaled  Burst  Depth  (m-kg~^) 


Figure  4.  Average  Scaled  Crater  Depth  Versus  Scaled  Burst  Depth 
for  Different  Soil  Characterizations  (Ref.  5)  Using 
Coefficients  From  Tables  I  and  11  (Ref.  6) 
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Figure  5.  Scaled  Crater  Volume  (Product  of  Crater  Radius  Squared  and  Crater 
Depth  Divided  by  Explosive  Charge  Weight  —  r2dc/W)  With  Standard  Deviation 
of  Measurements  for  the  29  Explosion  Events  in  the  DIRT  I  Test  (Ref.  7). 

The  bare-charge  events  (all  surface  detonations)  resulted  in  scaled  crater 
volumes  corresponding  to  the  dry/dry-to-moist  sandy  soil  categories  of 
Ref.  5.  Tube-delivered  projectiles  and  projectiles  statically  detonated  on 
the  surface  resulted  in  scaled  crater  volumes  which  were,  respectively,  401 
and  20%  smaller  than  these  resulting  from  bare  charges.  Buried  projectiles 
produced  larger  scaled  crater  volumes. 


Our  DIRTRAN  studies  show  that ,  st  nominally  fixed  soil  characteristics, 
a  vegetative  cover  produces  an  incremental  deposition  of  crater  material  at  the 
lip.  This  effect,  visibly  evident  in  Fig.  6,  can  account  for  the  different 
transmission  characteristics  exhibited  in  Figs.  1  and  2.  Vegetative  cover 
effectively  reduces  the  (elliptical)  crater  volume  (Fig.  7)  according  to  the 
equation 


V 

c 


r 


"d1 


where 

V  «  crater  volume  which  contains  only  soil  (i.e.,  no  sod) 
c 

r^.d^  ■  crater  radius,  depth  as  determined  in  accord  with  Refs.  5  and  6 
dg  “  sod  depth 


An  additional  correction  to  crater  volume  arises  in  connection  with  the 
explosion  of  bare  charges  vis-a-vis  casing  -  enclosed  charges  (e.g.,  artillery 
shells).  For  fixed  weights  of  explosive  and  terrains,  the  crater  size  directly 
depends  on  the  fraction  of  released  energy  which  is  coupled  into  the  ground. 

Not  only  does  this  fraction  vary  with  the  explosion  depth  as  noted  above,  but 
also  with  either  the  presence  or  the  absence  of  a  charge  casing.  For  155  mm 
artillery  projectiles  carrying  6.8  kg  of  TNT,  the  crater  volume  yield  (Fig. 5)  is 
approximately  60%  of  that  of  an  equal  weight  bare  charge  (Ref.  7).  This  evidence 
is  rationalized  upon  inspection  of  the  energy  balance  for  an  explosion. 

First,  let  us  consider  the  influence  of  tube  launch.  For  155  mm  projec¬ 
tiles,  the  kinetic  energy  associated  with  a  400  m/s  impact  velocity  is  14%  of 
the  explosion  energy;  its  contribution  to  increased  crater  volume  is  thus 
hardly  detectable  within  the  accuracy  of  the  crater  volume  measurements 
(+  20%).  Second,  let  us  consider  the  effect  of  casing  (Fig.  8).  Upon  explo¬ 
sion,  the  energy  is  partitioned  between  air,  soil,  and  (in  the  case  of  the 
projectile)  shell  fragments.  For  a  bare  charge,  the  partitioning  is  largely 
determined  by  the  compressibility  (elastic  modulus)  of  the  air  and  soil  media 
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^ *— I H  ^ 


Figure  6.  Graf  II  Test  155  mm  HE  Artillery  Round  Craters.  Depth  of 
vegetative  sod  cover  is  approximately  25  cm,  comprising  a 
large  portion  of  the  crater  volume. 
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Ft.  Knox,  Ky. 
Test  Range 


Grafenwohr,  Germany 
Test  Range 


Figure  7.  Effect  of  Sod  (Vegetative  Cover)  on  Crater  Volume.  The  sod 
depth  is  drawn  here  roughly  to  scale  for  the  Craf  11  test 
for  the  155  mm  HE  artillery  round  craters.  The  volume  of 
soil  (which  is  available  for  production  of  entrained  dust) 
is  seen  to  be  only  a  small  fraction  of  the  total  crater 
volume . 
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(4) 

shock  waves  in  the  two  media.  Detailed  cratering  calculations  show  that 
typically  5  -  10%  of  the  energy  is  coupled  into  the  6oil?  Therefore,  for  a 
bare  charge  exploded  at  the  surface,  90-95 %  of  the  explosive  energy  is  imme- 
diately  coupled  into  the  air.  Ultimately,  the  remaining  5-10%,  initially 
coupled  into  kinetic  energy  of  the  crater  ejecta,  is  also  transferred  to  the 
air  by  aerodynamic  arrest  of  the  particles,  as  depicted  in  Fig.  8a. 

The  explosion  of  a  projectile  leads  to  a  different  initial  partitioning 
of  energy  (Fig  8b);  the  following  arguments  indicate  that  an  appropriate 
estimate  for  the  explosive  energy  transferred  to  shell  fragment  kinetic 
energy  is  80%  with  20%  directly  transferred  to  the  air  thermal  energy.  Shell 
fragments  generated  by  the  upper  portion  of  the  exploding  shell  escape  ballis- 
tically  from  the  domain  of  influence  of  the  explosion.  By  contrast,  the  shell 
fragments  generated  by  the  lower  half  readily  impact  the  ground  and  thereupon 
transfer  their  kinetic  energy  partly  to  the  air  and  partly  to  the  soil.  Thesc- 
transfers  are  roughly  in  the  same  proportion  found  in  a  bare  charge  explosion 
because  of  the  impedance  matching  requirement  referred  to  above  and  discussed 
in  Ref.  4.  As  a  result,  the  energy  coupled  to  the  soil  is  approximately  40% 
less  than  the  equivalent  bare  charge  explosion;  consequently,  the  crater  is 
also  reduced  by  40%,  consistent  with  the  results  reported  in  Ref.  7.  The 
energy  ultimately  coupled  into  the  air  is  also  reduced  by  about  40%,  again 
consistent  with  the  data,  e.g.,  Fig.  9  where  the  development  of  buoyant  clouds 
produced  by  equal  weight  bare  and  enclosed  (projectile)  charges  are  shown. 
Clearly,  the  energy  loss  due  to  shell  fragments  results  in  smaller  buoyancy 
available  for  cloud  rise. 


The  same  estimate  can  be  obtained  upon  calculation  of  the  kinetic  energy  of 
the  ejecta.  For  example,  a  1  kg  charge  of  TNT  (4.7  x  IQ*1  J  energy)  exploded 
on  the  surface  of  dry/drv-to-moist  sand  (e.g..  White  Sands  DIRT  I  soil)  pro¬ 
duces  a  0.11  m^  crater.  Using  an  estimated  average  initial  ejecta  velocity 
of  60  m/sec  (Ref.  10)  and  a  specific  gravity  of  1.6,  the  total  kinetic  energy 
of  the  ejecta  is  then  3.2  x  10^  joules,  i.e.,  only  0.067  of  the  energy- 
released  by  the  explosion. 

These  estimates  lead  to  dual  predictions,  namely,  crater  sire  as  well  as 
buoyant  cloud  size,  which,  as  shown  below,  compare  favorably  with,  measured 
data . 


CLOUD  HEIGHT  (METERS) 


Another  consideration  concerning  the  dust  source  term  relates  to  the  fact 
that  only  a  fraction  of  the  crater  material  is  injected  into  the  atmosphere, 
the  balance,  having  a  small  ejection  velocity,  is  deposited  on  the  crater  lip. 
As  mentioned  above,  the  explosive  energy  coupled  to  the  ground  is  typically 
fixed  to  the  range  5-10%,  i.e.,  it  varies  at  best  by  a  factor  of  2.  How¬ 
ever,  the  crater  volume  for  different  soil  types  can  vary  by  a  factor  of 
40  for  surface  explosions,  as  evidenced  by  Figs.  3  and  4.  Thus,  by  conserva¬ 
tion  of  energy  for  a  given  charge,  the  kinetic  energy  per  unit  mass  of  the 
ejecta  associated  with  the  largest  craters  must  be  substantively  smaller  (by 
a  factor  of  the  order  of  20-40)  than  the  kinetic  energy  per  unit  mass  of  the 
ejecta  from  the  smallest  craters.  Correspondingly,  the  initial  velocities  of 
the  large-crater  ejecta  must  be  smaller  by  a  factor  of  the  order  of  4-6,  and 
the  fraction  of  the  ejecta  deposited  around  the  crater  rim  must  thus  be  larger 
than  for  small  crater  cases.  For  example,  for  explosions  in  granite  as  shown 
in  Fig.  10  (Ref.  10),  n0%  of  the  crater  ejecta  is  deposited  on  the  crater  rim 
and  40%  becomes  airborne.  From  the  kinetic  energy  considerations  stated 
above,  and  under  the  assumption  of  identical  impedance  matching,  the  airborne 
material  yielded  by  an  explosion  in  wet  cohesive  soils  could  be  reduced  by  a 
factor  of  4-6,  i.e.,  only  7-10%  of  the  crater  mass  becomes  airborne  while 
90-93%  is  deposited  on  the  rim.  These  percentages  only  represent  crude 
estimates.  However,  they  suggest  almost  an  order  of  magnitude  impact  on  the 
amount  of  material  becoming  airborne.  Thus,  either  systematic  cratering  cal¬ 
culations  or  systematic  experiments  focused  on  the  determination  of  crater 
rim  ejecta  are  essential  to  the  realistic  determination  of  soil  entrainment 
in  the  dust  clouds. 

2.1.2  Comminution/Agglomeration  of  Ejecta 

Even  if  the  factors  listed  above  are  taken  into  account  and  the  mass  cl 
dust  injected  into  the  atmosphere  determined  accordingly,  the  aerosol  is  not 
characterized  either  aerodynamically  or  optically.  Thus,  part  two  of  the 
present  problem  is  posed,  namely:  the  estimation  of  the  dust  particle  size 
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distribution  function.  On  aerodynamic  grounds,  this  is  important  because  all 
particles  having  diameter  D  1  100  pm  escape  ballistically  (cf .  Appendix  A) 
and  thus  never  partake  in  the  formation  and  subsequent  buoyancy/wind-controlled 
diffusion  of  the  obscuring  cloud.  On  optical  grounds,  the  determination  of  the 
distribution  function  controls  the  transmission  along  any  line  of  sight  through 
the  cloud  even  though  the  total  mass  and  dielectric  properties  of  the  inter¬ 
vening  aerosol  may  be  known.  Unfortunately,  the  "process  of  comminution  o: 
crater  material  has  not  been  described  either  theoretically  or 
experimentally"^' ^ 

Up  to  the  present  time,  dry-soil  sieve  analyses  have  been  used  as  a  mea¬ 
sure  of  the  distribution  function  to  be  encountered  in  explosion-generated 
clouds.  However,  prime  taeie  evidence  suggests  that  this  may  be  a  poor  assump¬ 
tion.  Specifically,  photographs  (Fig.  11)  of  the  ejecta  deposited  on  the 
crater  rim  indicate  the  presence  of  an  appreciable  amount  of  large  chunks  c : 
soil  ( a:  5-8  cm  diameter).  This  is  not  reflected  in  the  sieve  analyses,  su 
as  shown  in  Figs.  12  and  13  (which  represent  extreme  values  in  large  particles 
mass  fraction  indicated  by  the  sieve  analysis).  In  either  case,  no  more  •ha¬ 
ll  (at  most)  of  the  mass  is  associated  with  particles  having  diameter  I  _  ,  m . 
ir.  contradiction  with  the  photographic  data.  The  studies  of  optical  tra:.-: 
sior.  data  reported  in  Section  V  also  suggest  a  marked  depletion  of  I'.-.l 
particles  relative  to  the  amount  indicated  by  sieve  analysis.  In  fact,  :  . 
studies  indicate  that  particle  agglomeration  constitutes  the  largest  c.omu:  : 
xactor  in  the  determination  of  the  optical  effects  of  explosion-generat  c  d  c.:  : 
cloud.-;.  The  issue  warrants  immediate  attention.  Dynamic  soil  tests  which 
involve  the  failure  of  material  subjected  to  strong  wave  pulses  and  v.ich  cu:; 
provide  systematic  measurements  of  the  attendant  soil/sod  fragment  distribut  :  t 
function  are  recommended  as  indispensable  prerequisites  to  the-  systematic 
determination  of  the  airborne  dust  source  term. 

2 . 2  Initial  Dust  Clouds 

Inspection  of  data,  such  as  in  Fig.  14,  readily  reveals  that  the  initial 
dust  cloud  is  partitioned  into  two  parts,  one  buoyant,  one  nonbuovant .  Tin 


Effect  of  Dynamic  (Explosive)  Loading  on  Comminution/ 
Agglomeration  of  Soil  in  the  DIRT  I  Test  Series.  Not 
presence  of  large  ( « 5-8  cm  diameter)  crater  ejecta. 


ROT 


Figure  12.  Particle  Size  Distribution  Obtained  by  Dry-Soil  Sieve 
Analysis  for  DIRT  I  Test  Site.  Note  maximum  particle 
size  indicated  is  0.5  cm  (Ref.  8). 


29 


13.  Particle  Size  Distribution  Obtained  by  Dry-Soil  Sieve 
Analysis  for  DIRT  I  Test  Site.  For  this  sample, 
maximum  particle  size  indicated  is  1.8  cm  (Ref.  8). 
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Figure  14(b).  Example  of  Buoyant  Rise  Effects  on  the  Dispersion  of 
Explosion-Produced  Dust  Clouds.  The  effect  of  a  mild 
wind  which  causes  the  rising  cloud  to  follow  a  slanted 
path  is  evident.  Noticeable  too  is  the  persistence  of 
the  nonbuoyant  "dust  skirt"  which  is  near  the  ground 
and  which  has  significant  effects  on  ground-based 
EO  Systems  (Source:  DIRT  I  Test  Series). 
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geometry  of  these  clouds  can  In  principle  be  predicted  by  coupled  simultaneous 

time-dependent,  two-dimensional  calculations  of  the  soil  cratering  process 

and  the  aerodynamic  blast  propagation  process  on  time  scales  of  the  order  of 
-2 

10  sec  for  typical  artillery  explosions.  Such  calculations  are  not  avail¬ 
able  because  the  requirement  for  them  has  not  been  evidenced  heretofore.  One 
must,  at  present,  resort  to  data  analysis  and  correlations  guided  by  basic 
scaling  criteria  with  length  scales  proportional  to  the  cube  root  of  the 
explosion  energy  release. 


2.2.1  Buoyant  Dust  Cloud 

Whereas  the  explosion  energy  balance  discussed  in  Section  2.1  above 
indicates  the  amount  of  energy  coupled  into  the  air,  the  attendant  size  of 
the  buoyant  dust  clouds  can  be  estimated  on  the  basis  of  well-established 
correlations  of  hydrodynamic  blast  calculations  in  uniform  media  of  Infinite 
extent  (Ref.  12).  Thus,  the  initial  buoyant  cloud  radius  is  given  by 

1/3  <»</»> 1/3 


where  R^  is  in  units  of  meters;  Is  the  effective  energy  coupled  to  the 

air  in  units  of  kilograms  of  TNT;  p  and  p  are  the  densities  of  air  at  sea 

o 

level  and  at  the  blast  site,  respectively.  The  excess  temperature  AT  of  the 
cloud  over  ambient  air  follows  from  conservation  of  energy; 


AT 

T 

a 


jeff 


V  pc  T 
o  a  p  a 


-1 


eff 


V  p  c  T 

o  a  p  a 


0.57 


where 

*  energy  coupled  into  air  ■  kW^^ 

k  ■  energy  equivalence  of  TNT  ■  1.1  x  10^  kcal/kg  (Ref.  13) 

4  3 

V  -  initial  cloud  vclume  ■  r  »  K 
o  Jo 

3 

p  ■  ambient  air  density  ■  1.25  kg/m 
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hear  capacity  of  air  -  0.24  kcal/kg 


c  ■* 

P 

1  »  ambient  temperature  (assumed  a  288°K) 

.  tut.  initial  Just  cloud  temperature  T  is  1.57  T  (i.e.,  I  ■  45iLK  fur 
1  -  1:  k ;  . 

" he  celcular ions  presented  in  npf .  12,  however,  do  not  predict  the  le.o- 
rhe  ir.i*  •  '  buoyant  cloud  a-  rVe  e-  i  of  the  aerodynamic  transfer.!  , 

•  as  .  •  -w.  •.  i‘:b  th*  do  ii  sc  ;•  i  -  c6  pressure/vr’  nc-  ' ty  equil  ibrr- 

w  .  s  h  the  a;:.'  ii  ••  •  9tmc  .  pv  t  .  ...»  Eer  ■:  nu  d  a'tvi,  <  j  -.  •.  '  3c  ~  ' 

'■tie  v.  t •:.<  t.  -  »«'  .-.j  s  i e  P  .  Inspect!  o:  the  DIP'.  I  -loud  l.eig.  : 

o 

UiiU  -'Ref.  14)  for  sate  charges  exploded  l.:  th-'  s_tiact  dome.  ~t:  ate-  that  t 1 

he .t»  .  c.  ti.u  cloud  ctwitoi.'  ( z  *•  h-R  ,  where  r;  is  the  height  of  the  top  of  the 

•scaler-  witl  R  and  reaches  the  value  5R  at  the  end  o:  t  .  aerc._vr.ar.i 
o  o 

ui  ansie.-.t  time  t  «  t/R  ar  0.15  R  (for  t  in  seconds,  R  in  meter.-  -  i:  • 
do  o  d  o 

charge  is  not  exploded  on  the  surface,  the-  asymptotic  value  o:  t  >  e  centre  iu 
eight  clianges .  however ,  lirtt  •  prevent  a  completely  satisfactory 

c  'i  re  la  tier. .  Fo:  example ,  it  appears  that  for  exploi.'cs  o:  tuve-d.  .  in  ■  c  d 
artillery  project!!  ,  this  centroid  displacement  is  reduced  i  ? <  : o  (i.e., 
z  ■  R  ) . 

O 

*  direct  informatir:.  is  available  with  regard  to  initial  mas*  n  - 

ir-  ■  buoyant  cloud.  Because  typical  optica,  t  tans.-.ission  lir.e-of-.  :  f  : 

are  1  i 1 :  ?d  to  heights  less  than  5c,  they  intersect  the  buoyant  clot  1  ■  •. 

-.-■I’  fraction  of  a  second  (less  than  0.1  second'',  at  wh'ch  i»t 
ar<  op. .  ;ue  for  the  wavelengths  and  charge  sir-  employed  in  t  he  'IF!  1  .. 
f'.raf  11  series.  Hence,  only  lower  bounds  for  the  mass  loading  fa:  ’  •  d  if. 

:.  ii.e-i  from  t  he  so  data.  Obst-r  ration  o:  cloud  dissipation  at  late 
(cl.  he:,  t  ion  V)  indicate  that  perhaps  50%  of  the  small  size  comminuted 
agglomerated  ejecta  are  contained  in  the  buoyant  cloud.  This  only  represent: 
a  rough  estimate  which  must  be  improved  in  the  future.  Reduced  scale  experi¬ 
ments  with  smaller  charge  sizes  could  greatly  resolve  this  issue.  These 

experiments  should  be  executed  in  the  broader  context  of  answering  the  question 

of  initial  airborne  ejecta  mass,  comminution/agglomeration,  and  velocity. 


2.2.2  Nonbuoyant  Dust  Skirt 


The  existence  of  nonbuoyant  cloud  portions  is  forcefully  demonstrated 
not  only  by  photographic  data  (as  shown  in  Fig.  14)  but  also  by  the  long  dura¬ 
tion  of  optical  degradation  along  linea  of  sight  near  the  ground.  The  trends 
in  the  optical  data  discussed  in  Section  V  show  that  these  "dust  skirts"  dif¬ 
fuse  under  the  action  of  the  wind.  The  initial  dimensions  of  the  skirt  must 
again  scale  together  with  the  buoyant-cloud  characteristic  dimension  Rq. 
Inspection  of  typical  cloud  contour  maps  such  as  shown  In  Figs.  15  and  16 

(Ref.  15)  suggest  that  the  skirt  of  an  individual  explosion  can  be  approxi- 

* 

mated  by  upright  cylinders  of  radius  2R  and  height  2R  . 


The  data  in  Figs.  15  and  16  are  for  an  array  of  three  simultaneous  explosions 
oriented  perpendicular  to  the  line-of-sight  end  separated  by  a  distance 
d  ■  15m.  Thus,  the  apparent  initial  cloud  width  for  three  charges  is 
2(d  +  2RQ).  Accordingly,  one  predicts  for  Test  B-2  and  B-7  (charge  sizes 
6.8  kg  and  55  kg  TNT,  respectively)  apparent  widths  of  A5m  and  60m  which 
compare  reasonably  well  with  the  observed  values  of  AAm  and  54m. 
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Figure  15.  Computer  Generated  Contout  Plots  of  Dust  Clouds  Formed  in 
Event  B-2  of  the  DIRT-I  Test  0.5  sec  After  Detonation. 

Note  presence  of  three  overlapping  clouds  due  to  the 
three  6.8  kg  charges.  As  reported,  the  width  is  44  meters. 
The  separation  between  charges  is  15m.  Each  individual 
cloud  is  approximately  14m  in  diameter  (Ref.  15). 
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Figure  16.  Computer  Generated  Contour  Plots  of  Dust  Clouds  Formed  in 
Event  B-7  of  the  DIRT-I  Test  0.5  sec  After  Detonation. 

Note  presence  of  three  overlapping  clouds  due  to  the 
three  55  kg  charges.  As  reported,  the  width  is  54  meters. 
The  separation  between  charges  is  15m.  Each  individual 
cloud  is  approximately  24m  in  diameter  (Ref.  15). 
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3.  ATMOSPHERIC  DIFFUSION 


The  buoyant  cloud  evolves  at  first  like  a  thermal.  When  the  rate  of  dif- 
lus: .n  Jut  to  buoyant  rise  becomes  equal  to  that  associated  with  the  turbulent 
wind  field,  the  buoyant  cloud  development  becomes  controlled  by  the  latter. 

Just  skirt  evolution  is  controlled  throughout  by  wind-generated  turbulence, 
Q.-.ui-tative  dust  cloud  patterns  as  calculated  by  the  D1RTRAN  model  are  ex¬ 
hibited  in  Figs.  17,  18,  and  19  (Ref.  11).  Details  of  the  model  statement 
are  described  below. 


Buoyant  Cloud  Ri>  e  and  Diffusion 


This  process  is  described  by  the  model  of  Ref.  16  for  a  thermal.  With 
the  notations  l,  A,  and  V  for  the  cloud  radius,  surface  area,  and  volume 
respectively;  u  for  the  rise  velocity  (z-direction) ,  assumed  uniform  throughout 
the  cloud;  v  lor  the  horizontal  cloud  velocity,  assumed  equal  to  the  local  wind 
velocity;  AC  for  the  cloud  excess  potential  temperature  relative  to  local  poten- 

tial  temperature  i  ;  the  conservation  equations  of  mass,  momentum,  and  energy 

a 

take  the  form 


dh 

dt 


au 


(3-1) 


a 


1(10.  c 

dt 


(3-2) 


(3-3'' 


where  a  is  an  empirical  entrainment  coefficient  equal  to  0.23  according  to 
numerous  experiments  (Ref.  17)  and  g  Is  the  gravitational  acceleration. 
Equations  (3-1)  to  (3-3)  are  integrated  subject  to  initial  conditions  at  t  ■  0 
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Figure  17.  Wind  controlled  diffusion  of  nonbuoyant  dust  skirt  with 
gravitational  settling.  Time  sequence  of  contours  of 
line-of-s ight  integrals  of  dust  concentration 
(i.e.,  C£  values)  were  computed  by  the  DIRTRAN  model. 
The  three  contours  represent ,  respectively ,  0.1,  0.05, 
and  0.025  fractions  of  the  initial  average  CC.  Note 
the  effect  of  wind  shear.  (Ref.  11) 
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Superposition  of  Wind  and  Buoyancy  Controlled  Diffusion  of  Saal 
Particles.  The  time-sequence  of  contours  of  line-of-sight 
integrals  of  dust  concentration  (i.e.,  Ci.  values)  were  computed 
by  the  DIRTRAN  model  for  0.1,  0.05,  and  0.025  fractions  of  the 
initial  average  line-of-sight  integral  of  particle  concentra¬ 
tion.  Note  the  strong  wind  shear' effect  with  altitude,  an 
important  phenomenon  for  near-ground  line-of-sight  obscuration 
effects.  (Ref.  11) 
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for  bare  charges  detonated  at  the  surface 

for  tube-delivered  155  mm  artillery  shells  with 
fuzing  as  used  in  the  DIRT  1  test 


Upon  integration  of  Eqs.  (3-1)  to  (3-3),  the  height  h(t)  at  the  top  of  the 
cloud,  the  centroid  height  z(t),  and  the  cloud  width  d(t)  are  calculated  as 


z  (t) 


2(0) 


(3-4) 


h(t)  *  z ( t )  +  b(t) 


(3-5) 


d ( t )  *  2  b(t) 


(3-6) 


The  vertical  distributions  of  potential  temperature  0  (z)  and  wind  velocity 

a 

v(z)  arc  described  .according  to  atmospheric  stability  condition  as  detailed 
in  Appendix  B. 

Predictions  of  this  model  are  compared  with  data  in  Section  5. 

3 . 2  Wind-Dominated  Cloud  Diffusion 

The  process  of  wind-controlled  diffusion,  taking  into  account  the  actual 
vertically  sheared  wind  profile,  is  described  by  the  differential  equation  of 
diffusion.  The  action  of  turbulence  is  represented  on  a  long-time  average 
basis  by  the  gradient  transfer  relation.  The  dif fusivities  in  the  three 
* 

This  initial  velocity  can  be  extracted  from  the  point  source  solution  (Ref. 16) 
and  is  found  to  agree  with  the  data  for  cloud  rise  velocity  from  the 
DIRT  I  test. 


coordinate  directions  ere  prescribed  consistent  with  the  stability  character¬ 
istics  of  the  atmospheric  boundary  layer.  Aerosols  are  subjected  to  gravi¬ 
tational  settling  (at  their  sire-dependent  terminal  velocity)  in  addition  to 


diffusion. 


The  Initial  cloud  is  divided  into  horizontal  slices  of  circular  cross 


section.  Since  the  equations  are  linear  in  the  particle  concentration,  the 
solutions  for  each  slice  are  superimposed  to  produce  the  final  solution. 

The  individual  solutions  are  obtained  by  transform  methods.  Specifically, 
we  consider  moments  of  the  particle  concentration  function  C(x,y,z,t),  viz.. 


■as 

JJ x^yq  C  dxdy 


(3-7) 


where  x  is  the  coordinate  in  the  windward  direction.  If  the  diffusion  equation 


+  7 

3t 


(3-8) 


where  T  denotes  the  particle  flux  vector  components 


uC  -  D  , 

X  3x  ' 


(3-9a) 


D  f 

y  3y 


(3— 9b ) 


-  v  C  -  D  f 
g  z  3z 


(3-9c) 


(with  u  the  z-dependent  wind  velocity;  D  ,  D  ,  and  D  the  z-dependent  diffus- 

x  y  z 

ivities  along  x,  y,  and  z  directions;  and  v  the  particle  terminal  settling 


velocity),  is  multiplied  by  appropriate  powers  of  x  and  y  and  integrated  over 
all  x  and  y,  the  moment  equations  are  obtained: 


(~  -  ■q—  D  -  v  T  ■  pu  T  .  +  p(p-l)D  T  . 

\3t  3z  z  3z  g  3z/  pq  v  p— 1 , q  v  v  x  p-2,< 


+  q(q-l)D  T  _  =  R 

H  H  y  p.q-2  pq 


(3-10) 
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The  first  several  moments,  respectively,  describe:  the  mass  distribution  of 
dust  in  the  cloud  as  a  function  of  height  z  and  time  t  (p  *  q  -  0) ;  its  mean 
along-wind  (x)  position  (p  ■  1,  q  “  0) ;  its  mean  cross-wind  (y)  position 
(p  “  0,  q  «  1);  its  along-wind  spread  (p  ■  2,  q  ■  0) ;  and  its  cross-wind 

spread  (p  •  0,  q  ■  2")  . 

If  variables  are  made  dimensionless  with  respect  to  a  length  h  equal  to 
the  height  of  release,  a  velocity  [D^h/h] ,  and  a  time  £h2/Dz(h)J,  the  equation 
for  the  zeroth  moment  T^,  subject  to  a  diffusivity  «  Dz(h)zn,  initial 
cor.d  it  ions 


t  -  0,  T 


6(z-l) 


(3-11) 


and  the  boundary  conditions 


z  -  0,  z  (3T00/3z)  ■=  0 


( 3-12a) 


(3-12b) 


can  be  solved  in  closed  form  to  obtain  for  v  -  0 

g 


(l-n)/2 


h  (2-n) 


z2~n-H  2z 

T  (2-n)2t  1 (n-1)/ (2~n)  (J 

■  J 


l-n/2 


(2-n)  t 


(3-13) 


where  I,  denotes  the  modified  Bessel  function  of  the  first  kind  of 
(n-1)/ (2-n) . 

order  [(n-1 ) / (2-n)}  .  If  v^  i  0,  a  closed  form  expression  for  the  zeroth 
moment  may  be  obtained  by  heuristic  matching  of  two  asymptotic  solutions 
obtained  by  the  following  process.  The  equation 


3T00  n 

at  az  z 

m 


+  aT, 


(3-14) 


where  a  =  v  h/D  (h) ,  is  recast  in  terms  of  the  independent  variables 
g  z 


z  +  at  -  1 


to  obtain 


3T00  I  [  .  n  "M1 

ir  -  «  lK  '  “  '  ~J  ' 


which,  in  turn,  may  be  approximated  by 


8T00  S  f.  .n  8T00 

—  -  “J  "  ° 


in  C  >>  at ,  and  by 
3T 


00  n  3  T°°  .  o 

t—  -  (1  -  a)  - r-  -  0 

3T  3C2 


(3-15a) 


(3-15b) 


in  C  >>  ax.  The  solution  to  Eq.  (3-15a),  subject  to  the  initial  and  boundary 
conditions  (Eqs.  3-11,  3-12a,  3-12b) ,  is  readily  transposed  from  Eq.  (3-13), 
viz. , 


T<0)  /  r  T \  «  _ 

T00  ^  .3 


U+l) 


(l-n)/2 


h  (2-n) 


exp 


ii±n!2±i" 

(2-n)2t  . 


»  I 


(n-1) / (2 


_n)L  (2-n)2T  J 


(3-16a) 


The  solution  tq  Eq.  (3-15b),  subject  to  the  initial  condition  Eq.  (3-11),  the 
requirement  of  boundedness  as  (  ■*  +  ®,  and  the  restriction  (1  -  ax)  _>  0,  takes 
the  classical  Laplace  form 

T00^  ~-l/2  gxp  (3-16b> 

T 

where  ?  =  J  (1  -  aT)n  dx.  A  heuristic  composite  solution  to  Eq.  (3-14), 

0 

valid  in  z  >  0  and  x  ^  a  may  then  be  written  as 
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I00(C,T,a) 


(3-17) 


-  fU.CT)  T^U.T)  +  [l  -  f(C,OT)]  T^a.T) 

where  f(£,a-r)  denotes  a  suitably  selected  weight  function  (see  below)  satis¬ 
fying  the  requirements  f  -*•  0  in  £>>aT,  f  ■+•  1  in  £<«xt. 

A  judicious  choice  of  the  function  f(£,at)  is  obtained  upon  a  6tudy  of 
the  exact  solution  of  Eq.  (3-14)  for  the  special  case  n  «  1,  i.e.,  the  case 
of  neutral  stratification,  where  the  diffusivity  varies  linearly  with  the 
height  z.  In  that  case,  the  change  of  dependent  variable 

*oo(z,t)  ■  z  Too(z,t^ 

leads  to  the  equation 


-a  3*00  _  _3_  I"  1-a  3?0o1 

3t  3z  L  3z  J 


with  solution,  subject  to  the  initial  and  boundary  conditions  (3-ll)and  (3-12a,b), 


Too(z 


,t)  -  *aT00(*,t)  -  “P  [*  I1]  !«  V  '  °"18) 


Whereas  the  solutions  given  by  Eqs.  (3-16a,b)  are  also  known,  the  weight  func- 
t  or.  f(^,ai)  is  readily  determined  upon  comparison  of  Eqs.  (3-17)  and  (3-18). 

If  this  result,  strictly  obtained  for  n  •  1,  is  adopted  throughout  the  practical 
range  of  values  0  <  n  ^  4/3  associated  with  stable  to  unstable  stratifications, 
the  closed  form  solution  Eq.  (3-17)  for  the  zeroth  moment  becomes  completely 
defined.  The  availability  of  this  solution  substantively  simplifies  the  deter¬ 
mination  of  the  higher  moments  along  the  lines  presented  below. 

Returning  to  Eq.  (3-10),  we  note  that  the  same  differential  operator 

appears  in  the  equations  for  all  the  moments.  However,  while  the  equation 

for  Tqo  is  homogeneous,  those  for  the  higher  moments  include  known  forcing 

functions  R  .  Whereas  the  initial  conditions  are:  t  ■  0,  T_n  «  6(z-  1), 
pq  00 


and  T  •  0;  the  solution  Eq.  (3-17)  for  Tnr<  provides  the  Green's  function 
pq  uu 

required  to  construct  solutions  for  the  higher  moments  by  superposition. 
Specif icelly,  if  we  let 

Gtt.r.i.V.cO  =  f  o(t-t')][?-?']1/2  exp  - 

4(r-T*>  J 

*  (i  -  fft-c.  u+»2?? 

1  J  (2-n)  (t-t')  I  (2-n;  (i-t') 


(n-1) / (2-n) 


r  ,  l-n/2  *J 

2<LW)\n/2U/^i3 _ 

(2-n)2  (t-t')  J 


then,  the  solution  for  the  general  moment  becomes 


T  ^.‘O 
p.q 


(3-19) 


f  dC  f  R  (CV)  ca.C'.T.T'.o)  dx’  ,  (3-20) 

-1  i 

at 


vhich  can  be  readily  evaluated  by  numerical  integration. 

A  matrix  of  solutions  (3-20)  covering  the  ranges  of  the  parameters  n  and  a 
associated  with  Pasquill  stability  categories  A  through  F  (see  Appendix  B)  and 
dust  particle  sizes  D  100  urn  have  been  evaluated.  The  results  in  tabular 
form  are  included  in  the  DIRTRAN-I  code  together  with  interpolation  routines 
which  permit  determination  of  the  solutions  for  intermediate  values  of  the 
parameters  (Ref.  1).  By  this  process,  the  calculation  of  the  dust  concentration 
within  a  given  dust  size  bin  at  a  prescribed  location  and  time  achieves  the 
same  simplicity  as  the  evaluation  of  a  Gaussian.  Predictions  of  this  model 
are  compared  with  data  in  Section  5. 
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4.  OPTICAL  PROPERTIES  OF  DUST  CLOUDS 


The  DIRTRAN  model  calculates  the  transmittance  through  a  dust  cloud 

t 

according  to  the  relation 


T 


(4-1) 


where  l  Is  the  thickness  of  the  cloud  along  the  line  of  sight,  and  a  is  the 
mass  extinction  coefficient  which  depends  on  wavelength,  particle  size  dis¬ 
tribution,  and  material.  The  diffusion  calculations  described  in  Section  3 
yield  the  integral 


I 


-  Cl 


(4-2) 


where  C  is  the  average  concentration  along  the  line  of  sight.  The  model 
approximates  Eq.  (4-1)  by  the  relation 

T  *=  exp  (-  a  Cl)  (4-3) 

where  cT  is  estimated  from  laboratory  data  as  shown  in  Fig.  20  (Ref.  18)  and 
Inferred  from  ratios  of  transmittances  in  field  Experiments  (Refs.  7,  19,  20). 
In  principle,  if  the  chemical  cf  -'sition  and  associated  refractive  index  of 
the  dust  material  were  given  together  with  the  dust  particle  distribution 
function,  Kie  theory  could  be  used  to  systematically  calculate  a  (ignoring 
the  problem  of  particle  nonsphericity).  In  fact,  these  calculations  have 
been  carried  out  parametrically  by  Jennings  et  al.  (Ref.  21).  The  salient 
results,  shown  in  Fig.  21,  indicate  that  the  extreme  variation  in  extinction 
coefficient  for  soil-like  materials  as  a  result  of  refractive  index  variations 
due  to  different  chemical  constituents  is: 
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SCREENING  TEST-EA  "DIRTM-lqtMo«on-Peord» 21Do77 
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Figure  21. 


Aerosol  Extinction  Coefficient  vs  Wavelength  for  Aerosol 
Mass  Loading  (1.5  x  10~2  g/m3)  Characterized  by  Bimodai 
Particle  Size  Distributions  With  Geometric  Radii  0.5  and 
15  urn  and  Standard  Deviation  2.25  and  1.6  respectively. 

The  circles  indicate  extinction  values  for  typical  retrac¬ 
tive  indexes,  and  the  error  bars  indicate  the  extreme  ex¬ 
tinction  values  due  to  refractive  index  variations 
(Ref.  21  ). 
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(a)  negligible  (<5%)  for  wavelengths  X  4  1  pm 

(b)  +30%  for  X  w  4  pm 

(c)  +50%  for  X  «  10  yin. 

Jenning’s  calculations  for  the  variation  In  relative  extinction  coefficient 
versus  wavelength  are  also  supported  by  DIRT  I,  Graf  II,  and  Ft.  Sill  Dust/ 
Debris  measurements,  which,  within  +30%  scatter,  indicate  the  relative  values 
6hown  in  Table  4-1.  In  view  of  the  uncertainties  in  soil  comminution/agglom¬ 
eration  raised  in  Section  2  and  quantified  in  Section  5,  it  was  deemed  appro¬ 
priate  at  the  present  to  adopt  in  the  DIRTRAN  model  these  average  values  of  a. 


Table  4-1.  Relative  Mass  Extinction  Coefficients 


x  (urn) 

0.4-0. 7* 

0. 8-1.1 

3. 5-4.0 

8.5-12 

2100-3200** 

a(m2/g) 

1 

■SB 

mm 

0.5 

7  x  10-A 

* 

^visible 


550  m2/kg  (Ref.  18) 


These  values  are  for  Rayleigh  regime  particles  (relative  to 
mm  waves)  as  determined  from  the  asymptotic  attenuation  level 
(t  >  15  sec)  observed  in  the  DIRT  I  C-l  event. 


A  divergence  exists  between  the  extinction  coefficients  extracted  from 

Jennings  calculations , viz . ,a  ■  k/C  (where  k  is  given  by  Fig.  21  and  C  ■  1.5 
-2  3 

x  10  g/m  ),  and  those  reported  in  Fig.  20.  The  latter  values  are  approxi¬ 
mately  a  factor  5  higher.  Although  no  information  concerning  the  particle 
6ize  distribution  function  is  reported  in  Ref.  18,  the  Fig.  20  data  were 
adopted  in  the  DIRTRAN  model  because  they  represent  actual  measurements  of  a 
real  soil  sample. 
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5.  DATA  ANALYSIS,  MODEL  VALIDATION,  AND  DISCUSSION 


The  capability  of  the  buoyantly  rising  cloud  model  in  the  DIRTRAN  code 
is  iilustratedby  the  comparisons  shown  in  Figs.  22,  23,  and  24  for  data  from 
the  DIRT  I  test  series.  Test  B-2  (Fig.  22)  involved  a  3  x  1  array  of  6.8  kg 
bare  charges;  Test  E-3  (Fig.  23)  involved  a  3  x  4  array  of  6.8  kg  charges 
(155  mm  HE  projectiles);  and  Test  C-l  (Fig.  24)  involved  a  7  x  20  array  of 
6.8  kg  bare  charges.  The  atmospheric  conditions  in  the  three  cases  are 
respectively  consistent  with  Pasquill  categories  E,  C,  and  E.  The  departure 
between  theory  and  observations  for  C-l  at  80  seconds  after  detonation  is  due 
to  the  presence  of  a  temperature  inversion  (determined  by  acoustic  sounder, 
Ref.  22)  at  an  altitude  of  140m.  The  strong  sensitivity  of  the  predictions  to 
stability  category  is  also  shown  in  Fig.  23.  Thus,  the  accuracy  required  in 
near-ground  meteorological  measurements  must  be  such  as  to  provide  significant 
distinctions  between  vertical  temperature  and  wind  velocity  differences  asso¬ 
ciated  with  different  stability  regimes. 

The  capabilities  of  the  wind-controlled  diffusion  model  are  exhibited  in 
Figs.  25  and  26  for  events  B2  and  Cl.  The  wind  velocity  and  direction  (rela¬ 
tive  to  the  line  of  sight)  for  Event  B2  are  1.3  m/s  and  45°,  respectively; 
for  Event  C-l,  they  are  1.3  m/s  and  9°.  It  should  be  noted  that  these  wind 
velocities  represent  the  average  of  the  5-minute  average  readings  of  the  north 
and  south  meteorological  station  sites*  as  reported  in  Ref.  7.  The  wind 
directions  constitute  intermedin . o  values  between  the  5  minute  average  read¬ 
ings  at  the  two  sites,  readings  which  differ  by  approximately  100°  in  both 
events.  Because  calculations  show  that  the  time  for  initial  transmission 
recovery  is  controlled  by  wind  advection  of  the  dust  cloud  perpendicular  to 
the  line  of  sight,  the  importance  of  careful  meteorological  measurements  as 
near  as  possible  to  the  test  area  is  emphasized.  The  time  average  of  these 
measurements  should  also  be  made  on  time  scales  smaller  than  or,  at  worst, 

* 

These  sites  were  each  located  one  kilometer  from  the  test  area. 


CLOUD  HEI6HT  (METERS) 


COHPAR1SON  OF  DIRTRAN  NDDEl 
PREDICTIONS  KITH  CLOUD  HEI6HT  DATA 


(Source:  DIRT  I  Test,  October  1978) 


Time  After  Event  (Seconds) 


Figure  22,  Comparison  of  DIRTRAN  Model  Cloud  Height  Versus  Time 
Predictions  with  Data  From  Event  B-2  (three  6.8  kg 
bare  charges)  of  the  DIRT  I  Test  Series  . 
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Cloud  Height  (Meters) 


COMPARISON  OF  DIRTRAN  MODEL 
PREDICTIONS  WITH  CLOUD  HEIGHT  DATA 
(Source:  DIRT  I  Test,  October  1978) 


Time  After  Event  (Seconds) 

Figure  23.  Comparison  of  DIRTRAN  Model  Cloud  Height  Predictions  with 
DIRT  I  Data  for  Event  E-3  (twelve  statically  detonated 
155  mm  HE  artillery  projectiles).  Note  sensitivity  of 
calculations  of  cloud  height  to  atmospheric  stability 
category. 
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HEIGHT  (METERS) 


gM*"''-*-' 


COMPARISON  OF  DIRTRAN  MODEL 
PREDICTIONS  WITH  CLOUD  HEIGHT  DATA 

(Source:  DIRT  1  Test,  October  1978) 


Test  C-l  (7x20  Array  of  6.8  kg  TNT 
Bare  Charges) 

■  Data 

—“DIRTRAN  Calculation  ^ 


Temperature 
Inversion 

Layer  ^ 


'Wind  Speed  (2m  Height):  1.3m/sec 

Temperature  Gradient  (0  to  100m  Height):  -0.09°C 

Pasouill  Category:  E 


Time  After  Event  (Seconds) 

Figure  24.  Comparison  of  DIRTRAN’  Model  Cloud  Height  Predictions  with 
DIRT  I  Data  for  Event  C-I  (140  6.8  kg  bare  charges). 

The  existence  of  a  temperature  inversion  layer  at  14Qm 
(determined  by  acoustic  sounder,  Ref.  22)  results  in 
departure  of  model  predictions  from  actual  cloud  height. 


TRANSMITTANCE 


COMPARISON  OF  DIRT RAN  OPTICAL  TRANSMISSION 
MODEL  PREDICTIONS  WITH  DATA 


78  OCT  03 


EVENT  B-2 


0721  MOT 


ELAPSED  MINUTES 


Figure  25.  Comparison  of  DIRTRAN  Model  Optical  Transmission  Predictions 
with  DIRT  I  Data  for  Event  B-2.  Note  intermediate  recovery 
of  transmission  due  to  incomplete  merging  of  dust  skirts 
caused  by  the  individual  explosions  at  the  time  when  the 
midway  point  between  the  two  clouds  is  advected  across 
the  line  of  sight. 
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TRANSMITTANCE 


70  OCT  05  EVENT  C-1  0711  fIDT 


ELAPSED  MINUTES 


Figure  26.  Comparison  of  DIRTRAN  Model  Optical  Transmission  Predictions 
with  DIRT  I  Data  for  Event  C-1  for  Which  Transmission 
Recovery  Occurs  Only  After  Long  Times  After  Detonation 
(>10  Minutes) . 
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equal  to  the  obscuration  times.  Meteorological  measurements  should  also  be 
taken  at  different  sites  along  the  entire  length  of  the  line  of  sight.  The 
latter  requirement  is  manifested  by  the  optical  data  for  cases  B4 ,  B5,  B6, 

B7 ,  and  D1  where  the  "cloud  folded  back"  over  the  line  of  6ight  (Ref.  7). 

Additional  points  of  interest  in  Figs.  25  and  26  include  the  intermediate 
transmission  recovery  for  event  B2  (Fig.  25).  This  is  due  to  incomplete  merg¬ 
ing  of  the  dust  skirts  caused  by  the  individual  explosions  at  the  time  when 
the  midway  point  between  two  adjacent  clouds  was  advected  across  the  line  of 
sight.  The  existence  and  magnitude  of  this  intermediate  recovery  can  be 
reproduced  in  calculations  only  by  allowing  a  drastic  reduction  (quantified 
by  an  agglomeration  parameter  approximately  equal  to  seven)  in  the  concentra¬ 
tion  of  optically  active  particles  predicted  by  crater  volume  calculations 
and  dry-soil  sieve  analyses.  Such  reduction  can  be  explained  by  soil  agglomer¬ 
ation  under  dynamic  loading  (cf.  Section  2).  Two  additional  results  support 
this  view:  (1)  the  agreement  between  calculations  and  observations  for  the 
final  transmission  recovery  along  the  near-ground  line-of -sight  (Fig.  26); 

(2)  the  apparent  dissipation  of  the  buoyant  cloud  for  times  t  -  80  sec  after 
detonation  indicated  in  Fig.  22.  Whereas  only  one  comminution/agglomeration 
parameter  has  been  introduced,  the  congruence  between  predicted  time-dependent 
transmission  across  two  distinct  clouds  suggests  that  the  effect  ir.  question 
is  a  real  one  and  the  attendant  estimate  justifiable.  Giver,  the  magnitude  of 
the  parameter,  this  constitutes  the  dominant  influence  upon  optical  properties 
of  dust  clouds  produced  by  explosions  in  a  given  -soil.  To  be  sure,  the 
agglomeration  effect  extracted  in  the  B2  event  also  provides  good  agreement 
between  calculations  and  observations  for  event  Cl  (Fig.  26).  Noteworthy  in 
the  latter  case  is  also  the  success  of  the  model  at  very  long  times,  i.e., 
the  prediction  of  transmission  recovery  approximately  ten  minutes  after 
detonation . 

The  large  agglomeration  effect  exhibited  in  the  DIRT  I  data  analysis 
cannot  be  blindly  extrapolated  to  other  soils.  It  can  only  be  conjectured 
that  its  influence  on  dust  optical  properties  will  remain  large  for  other 
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soils.  However,  quantification  of  the  effect  remains  dependent  on  systematic 
soil  testing  under  conditions  of  (explosive)  dynamic  loading.  It  is  worth 
repeating  that  an  interesting  approach  may  be  provided  by  small  scale  experi¬ 
ments  with  reduced  charge  sizes  which  yield  partially  transmitting  dust  clouds 
at  all  times  subsequent  to  detonation.  Optical  diagnostics  at  several  wave¬ 
lengths  and  lines  of  sight  could  be  applied  to  such  clouds  to  obtain  spatially 
and  temporally  resolved  data  throughout  their  history. 
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6.  SUMMARY /CONCLUSIONS 


A  model  (D1RTRAN)  descriptive  of  time/space  dependent  optical  properties 
of  explosion-generated  dust  clouds  has  been  developed.  This  model  includes 
three  components:  a  dust  source  term,  atmospheric  diffusion,  and  optical 
transmission.  Data  analysis  indicates  that,  for  a  given  soil,  the  greatest 
uncertainties  lie  in  the  determination  of  soil  comminution/agglomeration 
under  the  dynamic  loading  produced  by  explosions.  When  an  empirical  parameter 
descriptive  of  this  effect  is  employed,  the  model  predicts  cloud  optical 
properties  which  compare  favorably  with  experiments  under  a  variety  of  explo¬ 
sion  patterns,  charge  types,  and  atmospheric  conditions.  A  computer  code 
(D1RTRAN-I)  which  embodies  this  model  has  been  delivered  to  the  Army  Atmo¬ 
spheric  Sciences  Laboratory  for  inclusion  in  the  Electro-Optics  Systems 
Atmospheric  Effects  Library. 
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APPENDIX  A 

AERODYNAMIC  ARREST  OF  PARTICLES 


The  retardation  of  a  spherical  particle  of  radius  r  by  aerodynamic  drag 
is  described  by* 


m 


dV 

dt 


paV 


V 


(A-l) 


where 


m 


V 

C„ 


mass  of  particle  *  —  rr  p^ 
air  density 
particle  density 
particle  velocity 
drag  coefficient 

2 

cross-sectional  area  of  particle  *  rr 


For  Reynolds  number  Re  «  2rV/vfl  <_  50  (where  v  is  the  kinematic  viscosity  of 

air),  the  drag  coefficient  may  be  approximated  by  C  »  24/Re,  while,  for 
5  v 

50  <  Re  <  3  x  10  ,  Cp  *  0.47.  Eq. (A-l)  can  be  readily  integrated  to  yield  for 
the  constant-drag-coefficient  regime  (C^  ■  0.47)' 


P 


at 

r 


(A-2) 


where 

»  velocity  at  time  t 

■  initial  velocity  at  t  ■  0 

a  •  constant  ■  3p  C_/8p 

a  D  p 


*  ,  -3 

We  have  neglected  the  inertia  terms  of  the  fluid  since  p  /p  ~  10 

a  p 
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i 


and  for  the  Stoke' s  Law  drag  regime  (C^  ■  24/Re) 


.  Vf 

In  — 
i 


9_  pava  t 

2  p  2 

P  r 


-bt 

2 


(A- 3) 


where 

b  “  constant  *  9p  v  /2o 

a  a  p 

The  distances  traveled  by  the  particle  having  initial  velocity  V  during  the 
time  intervalt  is,  for  the  constant-drag-coefficient  regime  (3x  10  >  Re  >  50), 


/  aV '  t 

s  “  —  In  ( 1  +  - 

a  \  r 


(A-i) 


£‘"W, 


a 


(A-5) 


and,  for  R  <  50, 


2 

6  *  r  vi  f1  * exp  (-bt/r2)] 


(A-6) 


V 


(A-7) 


The  total  distance  traveled  and  total  time  t^.  for  a  particle  to  come 

to  rest  are  calculated  as  follows.  Particles  with  radius  r  and  initial 

velocity  such  that  3  x  10"*  >  Re  >  50  are  retarded  according  to  Eqs.  (A-2) 

and  (A-4)  until  Re  “  2rV,/v  *  50  or 

f  a 


V 
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25v 

_ a 

r 


(A-6) 


From  Eq.  (A-5),  one  then  computes  the  distance  traveled  in  this  phase  of 
constant-drag-coefficient  as 


(A-9) 
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(A-10) 


From  Eq.  (A-2),  the  tlae  to  travel'  this  distance  is 


%-t) 


After  this  point  in  time,  In  Eq.  (A-8)  becomes  the  initial  velocity  for  the 

Stoke"s  Law  regime  (Re  <  50) ,  in  which  the  distance  traveled  s  and  time  t 

“3  *  ^ 

to  reduce  the  velocity  to  V^e  can  be  calculated  from  Eqs.  (A-6)  and  (A-3) 

respectively,  viz.. 


25v  r 

a 


(A-l 1) 


and 


3r 

b 


(A-12) 


Figures  A-l  and  A-2  show  parametric  curves  of  the  total  distance  sT  and  time 
of  arrest  tT  given  by 

sT  -  s,  +  6  (A- 13) 


‘l  +  C2 


(A-1A) 


2  -3  3  3 

for  ■  0.15  cm  /sec,  p^  *  1.2  x  10  g/cm  ,  and  »  2.5  g/cm  .  The  range 

of  initial  velocities  have  been  chosen  in  accordance  with  Ref.  10. 


Aerodynamic  Arrest  Distance  (m) 


Initial  Particle  Velocity  (m/s)y 


Particle  Radius  (m) 


Figure  A-l.  Aerodynamic  Arrest  Distance  Versus  Particle 
Radius  for  Four  Initial  Velocities 


Aerodynamic  Arrest 


Figure  A-2.  Aerodynamic  Arrest  Time  Versus  Particle  Radius 
for  Initial  Velocities  of  10,  50,  100,  and 
500  m/s.  (Curves  are  essentially  identifical 
for  these  range  of  radii  and  initial  velocities.) 
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APPENDIX  B 

DESCRIPTION  OF  THE  ATMOSPHERIC  BOUNDARY  LAYER  STRUCTURE 


Attention  is  restricted  to  the  constant  stress  layer  where  the  wind 
velocity  field  remains  planar.  Following  classical  notation  (Ref.  22),  the 
Monin-Obukhov  length  L  is  defined  as 


0  u. 


kV* 


(B-l ) 


The  profiles  of  wind  velocity  0(2)  and  potential  temperature  0(z)  for  unstable 
stratification  and  2  <  •  2  1  are  described  by 

A 


where 
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For  unstable  stratification  and  2  >  -  2L,  the  profiles  are 


u(r) 

6(2)  -  6(0) 


uo  -  0.95  u*  |L|1/3  2_1/3  k-1 


T  -  0.95  [L| 1/3  2_1/3 


(B-5) 


(B-6) 
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where 


T  ■  T 


*  [••  (6 
•  [••  ffl 


^m(-2)  + 


VH  (-2)  +  0 


O.75J 
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(B-7) 
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For  neutral  stratification,  the  wind  velocity  profile  is  described  by 

u(z)  -  -r-  in  (B— 9) 


For  stable  stratification,  the  wind  velocity  and  potential  temperature 
profiles  are  described  by  Eq.  (B-2)  with 
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-  7t 


=  -  11  f 


(B-10) 

(B-ll) 


for  z  <_  1 .5L,  and  by 

u  (z)  -  uQ  + 

0(z)  -  G(0) 

with 
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To  '  T*  [*-  (4f)  -  th  h-5) 

for  z  >  1.5L. 
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(B-l 3) 
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I 
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In  all  cases,  the  momentum  and  temperature  dif fusivlties  In  the  z  direc¬ 
tion  are  calculated  from  the  definitions 


2 

u* 

^zM  (du/dz) 


(B-15) 


ku*T* 

Dz0  “  (d0/dz) 


(B-16) 


with  du/dz  and  d0/dz  determined  consistent  with  the  profiles  given  above. 

In  accordance  with  Ref.  24,  the  corresponding  dif fusivlties  in  the  horizontal 

directions,  D  „  and  D  are  set  equal  to  5D  6D  w,  and  2D  ,,  for  neutral, 
xM  yM  n  zM  zM  zH 

stable,  and  unstable  stratifications  respectively. 

If  wind  and  temperature  measurements  are  given  at  two  heights,  the 
DIRTRAh'-I  code  inverts  the  wind  velocity  and  potential  temperature  profiles 
listed  above  to  calculate  case  by  case  the  quantities  u*,  T^,  z^,  and  L. 

Upon  this  determination,  the  wind  velocity  and  potential  temperature  at  any 
other  height  z  are  calculated  from  the  same  profiles.  If  wind  velocity  and 
temperature  are  measured  at  only  one  height,  the  DIRTRAN-I  code  requires 
that  Pasquill  stability  category  (A  through  F)  (Ref.  22)  and  roughness  height 
zq  be  prescribed  before  the  wind  velocity  and  potential  temperature  profiles 
can  be  determined.  With  each  stability  category,  there  is  associated  a 
specific  value  of  the  Monin-Obukhov  length  in  accord  with  Table  B-l  (Refs. 

22  and  25) . 

TABLE  B-l 

Monin-Obukhov  Length  vs  Pasquill  Stability  Category* 

- r 


Category 


L  (meters) 


